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Does the QCD plasma contain propagating gluons? 
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Comparison of two appropriately chosen screening masses of colour singlet operators in the pure 
glue QCD plasma indicates that at sufficiently high temperature it contains a weakly-interacting 
massive quasi-particle with the quantum numbers of the electric gluon. Still in the deconfined 
phase, but closer to T c , the same mass ratio is similar to that at zero temperature, indicating that 
the propagating modes are more glueball-like, albeit with a lower scale for the masses. We observe 
a continuity between these two regimes. 
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With the RHIC fully operational and busy taking data, 
certain questions about the treatment of the QCD plasma 
have become urgent. One of the most basic is about the 
modes of excitation in a plasma: does it contain weakly 
interacting quark and gluon quasi-particles, or some more 
complicated collective excitations? We show here that 
lattice computations yield detailed answers to this ques- 
tion. 

The Debye screening mass, m D , in pure gauge QCD, 
at temperature T > T c has an expansion in the strong 
coupling g of the form — 
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where the leading term g is well known, the second 
term has been extracted in perturbation theory fjj] and 
the non-perturbative coefficients b — 2.46 ± 0.15 and 
c = —0.49 ± 0.15 have been computed in a lattice sim- 
ulation of dimensionally reduced QCD ||. The QCD 
coupling g has to be evaluated at the scale 6.742T. Since 
T c /A— = 1.15 ± 0.05 @, g = 0(1) for T/T c ~ 3 or less, 
and the error due to the neglect of the g 4 term is about 
35<i%, where d is the coefficient of this term. As a re- 
sult, it becomes difficult to validate perturbation theory 
by comparing this screening mass to lattice data . 

In this work we address a prior question: at tempera- 
tures of interest to current and near-future experiments 
what mediates the longest correlations in the plasma? 
We examine this by comparing screening masses, m, ob- 
tained from correlations of two gauge invariant operators 
with different symmetry properties. They are chosen in 
such a way that one would be obtained by the exchange 
of two electric gluons while the other would need three, 
if indeed such gluons are the lightest excitations in the 
plasma. As a result, the two screening masses would be 
roughly in the ratio 3/2. 

Since screening masses involve the transfer matrix in 
a spatial direction, they are classified by the symme- 
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try group of the lattice sliced perpendicular to a spa- 
tial direction j|, || For the thermodynamics of a 
3+1 dimensional field theory realised on a hypercubic 
Euclidean lattice, since the Euclidean time direction is 
distinguished from the spatial directions, this is the group 
D4 x Z2(T) x ZiiCT), where D4 is the tetragonal group, 
Z 2 (C) is the charge-conjugation symmetry of the fields, 
and the Z 2 (T) factor arises from the symmetry t «-» — t 
II 0]- The transfer matrix can be block diagonalised in 
irreps of this group. 

Extensive thermodynamic quantities depend only on 
the lowest eigenvalue of the transfer matrix, and the 
phase structure is determined by the degeneracies and 
symmetries of the corresponding eigenvectors. In high 
temperature QCD, these are the scalar, irrep [0. 

In the confined phase, the unique ground state is scalar 
under the Z3 symmetry of the center of the colour SU (3) 
group. In the deconfined phase the ground state is three- 
fold degenerate and corresponds to the three irreps of Z3. 
On any finite lattice the degeneracy is lifted by an expo- 
nentially small quantity due to tunnelings between these 
states. In the sector, therefore, an unpbysical small 
"tunneling" mass, m T , may dominate the screening. We 
show later that we can control this and obtain the correct 
physical screening mass. 

Since the link variable, U ~ exp[J dxigA] ~ exp[m<7/l] 
(here a is the lattice spacing), one can turn this around 
and define a lattice gluon field of momentum k by the 
relation 

A,(k) = 1 e** [U»(x) - Ufa) - ImTr U,(x)] , (2) 

X 

where the sum is over all lattice sites in a slice, the com- 
ponents of k run over the set 2irl/N where 1 < I < N and 
N is Nt for the temporal momenta and N s for the spatial 
momenta (we assume an N t x N s 3 lattice). This defmit 1011 
gives an element of the SU (N) algebra which goes over 
into the continuum definition of the colour octet gluon 
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field in the limit of zero lattice spacing. The electric gluon 
A t -(0) is in the irrep, Since ® A~2~ = A^ + and 
A~2~ ® A~2~ <E> A~2~ = A~2~ , screening correlators in the 
colour singlet Af + and A^~ sectors would be dominated 
by two and three electric gluon exchange respectively. 
We discuss the influence of the magnetic sector later. 

In this work we use two classes of Af + operators — 
Wilson loops, specifically several linear combinations of 
the plaquette and the planar 6-link loop (sometimes 
called the fenster), and the trace of the real part of the 
Wilson line. For the A^~ operator we use the imaginary 
part of the trace of the Wilson line. The zero momentum 
projection is obtained as usual by summing over all sites 
in the slice. Since the Wilson lines are non-trivial un- 
der Zj, transformations, they are evaluated in the phase 
where the expectation value is real || . 

Due to the Z% symmetry of the vacua, there are only 
two tunneling masses of relevance — one in the vicinity of 
T c due to tunneling between the disordered state and any 
of the ordered states, and the other for all T > T c due 
to tunneling between any two of the ordered states (the 
latter are relevant only to operators which are non-trivial 
under Z3). Each mass has very specific dependence on 
the volume, V — a 3 N^N t , of a slice orthogonal to the 
direction of propagation — 



m r (\') = ( =L ) <'xp(-o-ri. 



(3) 



where C, a and a are constants. In a d-dimcnsional 
scalar theory a one-loop computation gives a — d/2 [ ^0[ . 
Tunneling arises when, on a finite system, simulations 
start exploring the non-Gaussian part of the free energy 
away from local minima. Clearly, an appropriate correla- 
tion function can give m T only when the order parameter 
distribution shows multiple peaks. A consequence of eq. 
(§) is that one can perform a finite size scaling study to 
check whether the lowest screening mass obtained is a 
tunneling mass. 

Details of our runs with the pure gauge Wilson SU (3) 
action are summarized in Table [j]. The critical coupling, 
(3 C , and its shift on finite lattices, is known for N t = 4 
with high precision Since the finite-size shift of 

(3 C on the smallest lattice, N s — 16, is less than 2 parts 
in 10 3 , the temperature scale is known with high preci- 
sion at T c . It is also known at similar precision at 1.5T C , 
2T C and 3T C from measurements with N t = 6, 8 and 
12. At other points the temperature scale is interpolated 
through the QCD beta function and has errors (|] , which 
are indicated in the table. Loop operators are measured 
at five levels of single link fuzzing Q. Cross correla- 
tions between all loops in the same irrep are measured 
and the lowest screening mass obtained by a variational 
procedure. Other details of measurements and analysis 
remain as in Q. 

At T — T c , the statistics collected are large enough 
that tunnelings between the deconfined and confined 
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24 


5.6500 


0.89 
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11370 


24 


5.6800 


0.97 
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12080 


16 


5.6908 


1.00 




20000 


20 


5.6918 


1.00 




20000 


24 


5.6920 


1.00 




20000 


16 


5.7010 


1.02 


(1) 


10000 


20 


5.7010 


1.02 


(1) 


10000 


24 


5.7010 


1.02 


(1) 


20020 


24 


5.7100 


1.04 


(1) 


10160 


24 


5.7200 


1.07 


(1) 


18010 


24 


5.8000 


1.27 


(2) 


32220 


24 


5.8941 


1.50 




15130 


24 


6.0625 


2.00 




31350 


24 


6.3500 


3.00 




31630 



TABLE I: For these runs on Nt = 4 lattices, measurements 
were taken every 5th sweep, where every sweep was done with 
a 3 hit pseudo heat-bath (except for the runs at T c , where 
measurements were taken every 10th sweep, each with a 6 hit 
pseudo heat-bath). In each case, an initial 1000-3000 sweeps 
were discarded for thermalization, and masses were extracted 
by a jack-knife analysis with 100 bins. 




FIG. 1: A~l + screening masses at T c and their dependence 
on the lattice size. The lowest mass (unfilled circles) is a 
tunneling mass, as evidenced by the good fit to the form in 
eq. Also shown are estimates of the second variational 

mass (filled circles) and a comparison with a measurement of 
the lowest mass at T/T c = 0.97 on a 4 x 24 3 lattice (horizontal 
band). Data for N s = 8 is from 



phases occur many times, as do those between different 
deconfined phases. As a result, it is only to be expected 
that the lowest screening mass that we can extract in the 
A^ + sector is the tunneling mass. Evidence for this is 
the good fit to eq. (J3J) shown in Figure 0; the fit gave 
X 2 — 1.005 per degree of freedom. 

To stabilise the pure phases we move away from T c . 
The distance from criticality, A/3 = \(3 — f3 c \, needed to 
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remain in a single phase for arbitrarily long runs depends 
on N s : A/3 can decrease exponentially with N s . In the 
confined phase we have estimated m(T c ) by a measure- 
ment with A/3 = 0.0125, corresponding to a 3 ± 1% shift 
of T below T c . A 4 x 24 3 lattice was observed to stay 
in the confined phase throughout the run. The second 
variational level was seen to correspond to the screening 
mass in this phase. Our measurement of the A^ + screen- 
ing mass at T/T c = 0.97 gives m{T)/T c = 3.41 ± 0.08 
(at T/T c = 0.89 we have m(T)/T c = 3.52 ± 0.07). 
In comparison, a zero temperature measurement of the 
scalar glueball mass at a similar lattice spacing gives 
m(T = 0)/T c = 3.93 ± 0.05 §J|. At T = Q, glueball 
masses depend strongly on the lattice spacing, whereas 
mass ratios are less sensitive. In this phase, we may then 
expect that the ratio m(T)/m(T — 0) may scale better 
than the ratio m{T)/T as the lattice spacing goes to zero. 
Our measurements show that m(T)/m(T = 0) ~ 0.8 in 
the scalar channel near T c , in agreement with recent mea- 
surements from temporal correlators fl4|| . 

Another method of extracting the physical correlation 
length is to pick out the configurations in which the 
whole lattice is in a single phase. We do this through 
the distribution of the action density, which has peaks 
corresponding to each of the phases. By restricting mea- 
surements of correlation functions to configurations with 
some S > S cu t, suitably chosen, we can isolate the de- 
confined phase and measure the physical screening mass. 
This procedure on the largest lattice gave m(T c )/T c 1 
in the deconfined phase. However with increasing S cu t 
statistics become poorer; consequently, errors increase 
rapidly, and it is hard to quote a more precise value. 

We worked in the deconfined phase by taking A/3 > 
0.0085 for N s = 16, 20 and 24, corresponding to mov- 
ing off from T c by 2 ± 1% or more. None of our T > T c 
runs showed any tunnelings between the different ordered 
states except the run at T/T c — 1.02 on the 4 x 16 3 lat- 
tice, where we found one tunneling event between two 
deconfined vacua (which forced us to use only the corre- 
lations of loop operators here). Within the precision of 
our measurement, the screening masses at T/T c = 1.02 
are independent of the lattice size: an indication that 
they are not tunneling masses. For all our measurements 
of the Af + , shown in Fig. 0, the masses extracted from 
the Wilson loops and those from the Wilson line agree at 
the 95% confidence level. 

The rough agreement between the measurement at 
1.02T C and that at T c in the deconfined phase, displayed 
in Figure |[ should also be noted. We also draw attention 
to the feature that the screening masses, when expressed 
in units of T, dip near T c . If the dip persists, then it 
cannot be understood in the context of the perturbation 
theory leading to eq. ([!]). An attempt to capture this 
effect in a model has been made in jl6) . We plan a more 
detailed study in both phases of the region near T c . 

In pi screening masses have been extracted from the 
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FIG. 2: Af + and screening masses as a function of tem- 
perature on 4 x 24 3 lattices. Note the discontinuity in the 
Af + screening mass at T c . 



free energy change due to the addition of a static quark 
pair in the plasma, i.e., the logarithm of the point- 
to-point Wilson-line correlation. Our measurements of 
m(T)/T for the scalar in the deconfined phase are com- 
pletely compatible with their results. This agreement 
is non-trivial since we extract the screening mass from 
zero-momentum screening correlators. 




FIG. 3: Ratio of the screening masses in the A 2 and A^ + 
sectors as a function of temperature. 

Our results for the screening mass are also shown 
in Figure || jl5| . In Figure || we have displayed the same 
data in the form of the ratio of the two screening masses. 
Clearly, already at temperatures a little above 1.25T C , 
the ratio is close to the perturbative value of 3/2, ap- 
proaching it from above. It has been shown earlier j(| 
that no mass assignment for magnetic gluons simultane- 
ously satisfies the and screening masses, both 
of which should arise due to the exchange of two such glu- 
ons. If these colour singlet channels are treated as collec- 
tive excitations, then their large screening mass prevents 
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them from contributing through a pair exchange to the 
correlators. Even if we choose to disregard this ar- 
gument, and insist on the existence of magnetic gluons, 
their screening mass is very high (in this temperature 
range) simply because of the large screening mass of the 
B± . The simple fact that the Af + has the smallest 
screening mass protects the argument of electric gluon 
exchange when the ratio m(v4 2 '~)/m(A^ + ) is seen to be 
3/2 as in Figure |[ 

We have not extended our measurement of the 
screening mass below T c since previous studies have 
shown that in the confined phase the degeneracies of the 
screening mass are those expected from the T = sym- 
metries of the transfer matrix ||. At T = the A^" 
is just one component of a very heavy vector glueball 
{T^~ in the usual T = notation), and the measure- 
ment of its screening mass would require large numbers 
of operators and a significantly larger data set. It is in- 
teresting, though, that T = measurements, made over 
a broad range of lattice spacings, show that the ratio 
of this vector mass to the scalar is roughly 3 fl7|] . The 
ratio m(Af + )/m(A2~) for T > T c seems to interpolate 
smoothly between the T = and the infinite temperature 
values, crossing over rapidly from one regime to another 
in the temperature range between T c and 1.25T C . 

It is interesting to recall that the entropy density of the 
plasma is relatively low for T/T c < 1.1, and begins to sat- 
urate only for T/T c > 1.25 Q. While thermodynamics 
is more sensitive to short distance modes and screening to 
long-distance modes in the plasma, the two observations 
can nevertheless be due to a single cause, if the partition 
function for T/T c < 1.1 were essentially saturated by the 
colour singlet modes that we have seen. At sufficiently 
high temperatures, the propagating electric gluons can 
contribute significantly to the entropy 

In summary, we have shown that the screening in the 
scalar (A^ + ) sector is discontinuous across the pure QCD 
transition at T c and falls abruptly by a factor of nearly 
3 in going from the confined to the deconfined phase. 
At larger temperatures, this screening mass rises faster 
than linearly in T. Above T c , the ratio of the A^ 
to the scalar screening mass interpolates smoothly from 
a value close to the T = ratio to the value 3/2 ex- 
pected when these correlations are saturated by electric 
gluon exchange. The agreement with the latter value 
for T > 1.25T C indicates the presence of weakly cou- 
pled, massive quasi-particles with the quantum numbers 
of the electric gluon as the lightest excitations of the QCD 



plasma. Closer to T c the mass ratio is more similar to its 
T = value. Taken together with the smallness of the 
entropy density, this indicates that excitations are more 
nearly glueball-like. 

We would like to thank Rob Pisarski for correspon- 
dence on the temperature variation of the A^ + screening 
mass. SD is supported by GSI under project number 
BI-KAR. 
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